Mesoscopic conductance fluctuations in a coupled quantum dot system 
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We study the transport properties of an Aharonov-Bohm ring containing two quantum dots. One 
of the dots has well-separated resonant levels, while the other is chaotic and is treated by random 
matrix theory. We find that the conductance through the ring is significantly affected by mesoscopic 
fluctuations. The Breit-Wigner resonant peak is changed to an antiresonance by increasing the ratio 
of the level broadening to the mean level spacing of the random dot. The asymmetric Fano form turns 
into a symmetric one and the resonant peak can be controlled by magnetic flux. The conductance 
distribution function clearly shows the influence of strong fluctuations. 
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The quantum dol*^ (QD) is an ideal system to study 
the phase coherence of quantum mechanical wave func- 
tions. Such effects can be explained by the interference 
of different pathways, induced by elastic scattering from 
an irregular boundary and/or impurities. When a dot is 
connected to leads, there is a strong overlap between the 
dot and leads and we must treat the whole as a quan- 
tum system. In fact, a recent numerical calculation of 
a chaotic dot shows that there are scars that connect 
between the leads £ 

The coexistence of a direct path and discrete levels in 
the dot induces a prominent effect, the Fano effect^ It 
was shown in QD systems^ that asymmetric Fano peaks 
can be controlled by gate voltages and magnetic fields. 
The work in Ref. demonstrated that an Aharonov- 
Bohm (AB) ring is a suitable system to study this effect. 
A dot is embedded in one of the arms, and the ring ge- 
ometry is utilized as an interferometer. This effect has 
also been observed in other experiments, including a mi- 
crowave cavitj*£ and optical absorption^ 

The situation changes drastically if we take into ac- 
count mesoscopic fluctuations. A random distribution 
of levels in the dot leads to sample-to-sample fluctu- 
ations of the conductance. In contrast to bulk systems, 
such fluctuations cannot be neglected in QD systems, and 
the conductance distribution has a broad non-Gaussian 
shaped It was observed in a QD system 10 and in a chaotic 
microwave cavityiii Clerk et al. developed a statisti- 
cal theory of Fano peaks assuming a random distribu- 
tion of peaks>i£ These authors focused on resonance- 
to-resonance fluctuations of the asymmetric Fano form, 
rather than conductance fluctuations. 

In this paper we develop a statistical theory of the AB 
ring with two QDs in the arms. In addition to a reso- 
nant dot in one of the arms, a random dot is connected 
to the other arm and is treated by random matrix theory 
(RMT)pi^ The orbit through the resonance is correlated 
to those through the random levels, and strong fluctua- 
tions can be observed in the conductance. We show that 



the Breit-Wigner and Fano resonant forms are no longer 
maintained in the averaged conductance. 

We focus on the following two issues. First, we ex- 
amine the properties of those orbits that contribute to 
the conductance. Motivated by the work of Ref. 0, we 
take into account the direct nonresonant path in the ran- 
dom dot as well as the Breit-Wigner resonant path in the 
regular dot. We show that these two channels give quali- 
tatively different contributions to the conductance in the 
presence of random levels. Second, the level broaden- 
ing of the random dot. For an open dot, the broadening 
can be larger than the mean level spacing due to strong 
coupling to the leadsi^ We systematically change the 
broadening from small to large values to elucidate how 
this parameter affects the results. 

Our system is defined by the internal Hamiltonian for 
the two QDs and their couplings to the left and right 
leads. It is depicted in the upper left inset in Fig. QJa). 
We assume that one of the dots (dot 1) has regular res- 
onant levels and the level spacing is much larger than 
the level broadening. In this case, each of the levels can 
be treated independently. The other dot (dot 2) is rel- 
atively large and has many levels distributed randomly. 
As is known from scattering theory (see, e.g., Refs. Ill2l 
and ll4|) . the S matrix of the system is written as 



5 = 1 — 2mw^ - 
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E-H 



(1) 



where H = H\ @ Hq, denotes the Hamiltonian for the 
QDs, w is the dot-lead coupling, and E is the (Fermi) 
energy. We adopt a single-level Hamiltonian Hi = E\ for 
dot 1. H% for dot 2 is a member of the Gaussian unitary 
ensemblei^ and its size N is taken to be infinity to find 
the universal limit. We assume that the S matrix is a 
2x2 matrix, which means the left and right leads have a 
single channel respectively. The conductance, measured 
in units of 2e 2 /h, is calculated from g = (|>Si2| 2 ). 

In the present RMT approach, there is no need to know 
the full form of w. The matrix w appears in the S matrix 



2 



formula as w^Gw, where G = (£'+ — Zf) 1 . After the 
averaging, the matrix structure of G is lost and 2x2 



matrices 7, 



7TW 



appear in averaged quantities. 



Here i = 1,2 label the dot and it^ 1 ) (k/ 2 ) ) is a 1x2 (A x2) 
matrix. Assuming a symmetric coupling with respect to 
the left and right leads, we arrive at the general form 



71 = T 



-lip 



AT, 



72 



<7C 



,(2) 



where 1^ is the level broadening of the dot i and ip is 
the AB fluxii^ The real parameter a (0 < a < 1) repre- 
sents the strength of the nonresonant direct path since 
the off-diagonal part of 7 contributes to the conductance 
directly^ 

The averaged S matrix is calculated from the relation 
(1-<S»/(1+(S» = (1-S 1 )/(1+S 1 )+(1-(S 2 »(1+<S 2 », 
where Sip are the S matrices of the dot 1,2 respec- 
tively. In RMT the magnitude of the Gaussian fluc- 
tuations determines the mean level spacing A. Then 
the averaged Green's function of the dot 2 is given by 
(G) = —iir/NA, where we take the limit N — > 00 
and neglect the real part of G. In this limit we ob- 
tain (5 2 > = (1 - 7r 72 /AA)/(l + 7T7 2 /AA) while Si = 
(E- Ei- vy\)/(E -Ei+vyi). (S) is determined by the 
following four parameters: a, ip, e = (E — E\)/Ti 1 and 
X = 71T2/ A. e measures the distance from the resonance 
point of the regular dot and X represents the ratio of the 
level broadening to the level spacing for the random dot. 
When X ^> 1, the level spectrum is continuous. The 
use of RMT implies our consideration is restricted to the 
energy scale much smaller than the Thouless energy. 

If we disregard quantum fluctuations, the conductance 
is given by go = |(S'i2)| 2 , which we call the principal part 
of g. It is given by the Fano form 



a 2 A 2 
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(l + a)X/2] z [l + (l 

\E-E 1 +qT 1 



a)X/2Y 
2 



(E-E!f 



[l + (l-acos2tp)X/2Y 



,(3) 



^p2 

[l+(l+a)X/2Hl + (l-a)X/2] 2i 1 

where q — ie~ 2%tp jaX is the Fano parameter. This pa- 
rameter is complex and becomes purely imaginary when 
the AB flux ip is zero, which is contrasted with the ex- 
periments for clean systems^ where a real q has been 
observed in the absence of a magnetic field. We also 
note that the result of Eq.JSJ holds regardless of the 
choice of the universality class because the resonance is 
not treated randomly, in contrast with the approach in 
Ref. 0. When a = 0, \q\ — > oo and the result reduces 
to the standard Breit-Wigner form. The presence of the 
random dot leads to a reduction of the level broadening 
and the conductance by a factor of 1/(1 + A/2) 2 . 

We now consider mesoscopic fluctuations of the con- 
ductance Sg = g — go. We calculate these using the 
method of supersymmetry^ which allows us to derive 
the nonlinear a model 

T/2 Aa + aA \ 
l-T/2 2 J ' { ' 



F = -str 8 In ( 1 



where the 4x4 supermatrix a parametrizes the saddle 
point manifold and satisfies a 2 — 1^ str denotes the 
supertrace and A = diag(l, — 1) in retarded-advanced 
space. The 2x2 matrix T defined by T = 1 - (5 , )(S , } t 
is called the transmission coefficients. This form of the a 
model is well known as a standard model of a single ran- 
dom dot. The only difference is in the T matrix. The fact 
that the a model is written in terms of T only demon- 
strates the universality of the correlation functions of the 
S matrix elements. For instance, Sg in Eq.Q is a func- 
tion of T. This is to be contrasted with the result of 
Eq.(J3J), where such T universality does not hold. In the 
present model, we find the eigenvalues of T at a = 



Ti 



2X 



(1 + X/2Y 



To 



2X 



(l + X/2) z + l/e 2 



(5) 



We see that T x is independent of e. 

We present the analytical results for the conductance 
together with numerical ones. Numerical calculations are 
performed using the random S matrix model^i where the 
internal structure of the S matrix is disregarded and ran- 
domness is imposed directly on S. For the distribution 
of S, we use the generalized circular unitary ensemble 
(CUE) based on the Poisson kernel^ 



P(S)dfi = 



dct(l- (5)(5)t) 2 
|det(l-S(S)t)| 4 



(6) 



where d/i is the measure of the CUEii^ This is the max- 
imally randomized distribution under the condition that 
the average value is (S) . Following the previous approach 
for regular systems^ we separate the system into the dots 
1 and 2, and the left and right forks that connect the dots 
to leads. Choosing the fork S matrix in a symmetric form 
we find the transmission through the ring expressed by 
the S matrix Si^ for each dot. S2 is treated statistically 
using a Poisson kernel. This random S matrix approach 
is equivalent to the random Hamiltonian approach if {S) 
is chosen properly and N — > 00,— We use the same ex- 
pressions for Si and (S%) as in the random Hamiltonian 
approach. 

First we consider the X dependence of the conductance 
at a = 0. When e — > 00, where the regular dot is detached 
from the random dot, we have T x = T 2 = 2X/{l + X/2) 2 , 
and recover the known result 20 



T T 2 
.90 = 0, Sg=^- + ^. 



(7) 



At the resonance point e = 0, T x = 2X/(1 + A/2) 2 , 
T 2 = 0, and we obtain 



1 



90 



(1 + X/2Y 



5g = - 



X/2 



(1 + X/2Y 



(8) 



In Fig. [IJa), X dependence of the conductance is shown 
for several values of e. go shows a peak at A = while 
Sg takes a maximum at A = 2, as shown by the thin 
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FIG. 1: (Color Online) Conductance vs X — 7TT2/A (a) and 
e = (E- Ei)/Ti (b) at a = 0. The thick (thin) lines are the 
total conductance g (the principal part go). Upper left inset 
in (a): Sketch of the sample. The parameters characterizing 
each dot are shown in the figure. Upper right inset in (a) and 
(b): Fluctuation part 8g = g — go- 



lines and the inset in Fig.QJa), respectively. As e — » 00 
go (Sg) is monotonically decreasing (increasing) and the 
result rapidly approaches Eq.JJJ). 

e dependence of the conductance is shown in Fig.[T^b). 
A resonance peak appears at e = 0, reflecting transport 
through the regular dot. This peak structure, however, 
changes qualitatively as a function of A. For small A the 
peak is convex and the peak height decreases on increas- 
ing X. When X = 2, g is independent of e. Increasing X 
further, we find that the peak turns into an antiresonancc 
and g decreases monotonically. The result for X = 2 cor- 
responds to that of the CUE because (52) = 0, and the 
Poisson kernel P{S2) becomes unity. 

The obtained results are for a model with a resonance 
at a = 0. Now the question is whether they are character- 
istic of the resonance model. For comparison, we discuss 
the direct reaction model with a^O and e — > 00. Both 
models are similar in that the conductance go becomes 
finite in the absence of fluctuations, go becomes maxi- 
mum at a = 1, which is similar to the situation at e = 
in the resonance model. In the direct reaction model, the 
eigenvalues of the T matrix are given by 



Tl,2 = 



2(l±a)X 



[1 + (1 ± a)X/2Y 



(9) 
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FIG. 2: (Color Online) Conductance vs e at a = 0.7 and 
ip = 7r/8. The thick (thin) lines are the total conductance g 
(principal part go). Inset: Fluctuation part Sg. 



T2 goes to zero as a — > I as in the limit e — > in Eq.JSJ. 
In this limit, we find T x = 4X/(1 + A) 2 , T 2 = 0, and 



x- 



9o 



(i + xy 



Sg 



Ti 
4 



A 



(i + x) 



(10) 



Note that T\ dependence of 5g is the same as in Eq.ljHJ 
but the A dependence is different. We have found numer- 
ically that the peak is maintained for an arbitrary value 
of X and that there is no antiresonance, in contrast to 
the resonance model. 

As we have shown in Eq.@, when we consider both 
the resonance and direct reaction channels, the Fano ef- 
fect appears in g . In Fig. [21 a typical numerical result 
of the conductance is shown at a = 0.7 and <p = n/8. 
We observe that an asymmetric form is obtained for the 
fluctuation part as well as for the principal part. How- 
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FIG. 3: (Color Online) Conductance distribution at a — 0. 
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ever, these asymmetric peaks work to compensate each 
other, resulting in a symmetric g that has a form sim- 
ilar to the previous result at a = 0. The difference is 
that the conductance as a whole is enhanced due to the 
direct path contribution. Far from the resonance point, 
the conductance is independent of (p, while it is sensitive 
at the resonance point. The resonance (antiresonance) 
is amplified at (p = (7r/4) where the Fano parameter q 
is pure imaginary (real). When q takes a complex value 
(when < ip < 7r/4), a dip in the resonance peak is 
formed at the intermediate values of X (See the plot of 
X = 2 in Fig. |2J). 

As we have shown, Sg are of the order of go and 
the mere calculation of the averaged conductance is not 
enough to characterize the system. We show the nu- 
merical results for the conductance distribution function 
P(g) = (S(g-\S 12 \ 2 )) in Fig.|3| For e > 1 we find the re- 
sults in Refs. IqI andll7l expressed in terms of T\ in Eq.JSJ). 
For finite value of e, the T universality does not hold 
and the results essentially depend on X and e. When 
X < 2, the peak representing the resonant conductance 
moves from g — to 1 as e decreases. When X > 2, the 
resonant conductance is suppressed and the distribution 
function shows the strong influence of the chaotic scatter- 
ing. When e = 0, two peaks appear at g — and 1. The 
peak at g — is larger (smaller) than that at g — 1 when 
X > 2 (X < 2), which clearly shows the coexistence of 
the contributions from both regular and chaotic dots. At 
X = 2, P(g) is always symmetric and is consistent with 
the results of the averaged conductance. The curve for 
e = is well fitted by the function P(g) — \/it\Jg(\ — g). 



We remark on the effect of dephasing. ft can be simply 
studied by introducing an imaginary part to the energy 
E — > E + i/2r. The substitution of this to Eq.© leads to 
the reduction of the conductance. The effect on the fluc- 
tuation part is to add F T = (Ar) _1 str oA to the a model. 
Numerically we have observed that the fluctuation part 
is strongly suppressed by this effect while the principal 
part shows a small reduction. This means that the reso- 
nance is preserved at any X. We also anticipate that an 
asymmetric Fano resonance can be observed, since the 
cancellation of the asymmetry between the principal and 
fluctuation parts becomes incomplete. The quantitative 
estimate of the dephasing effects using the method in 
Ref. is necessary to compare with the experimental 
results, as was done in Refs. Ilfl andllll for the single-dot 
systems. A detailed study will be reported elsewhere. 

Another interesting problem is the ensemble depen- 
dence of the result. Our numerical calculations using 
the orthogonal and symplectic ensembles show that the 
distribution at the resonant point (the lowest graph in 
Fig. |2J is independent of the ensemble, ft implies the 
distribution is determined by some universal mechanism. 

In conclusion, we have developed a statistical theory 
for an AB ring system with regular and chaotic QDs. The 
conductance and its distribution are strongly influenced 
by the mesoscopic fluctuations of the chaotic dot and the 
position of the resonance peaks. 
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Society for the Promotion of Science. 
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